Stability analysis of kinked DNA in T{K, r) model 
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Abstract 

We phenomenologically analyze short DNA rings' stability by discussing the second variation 
of its elastic free energy. Through expanding the perturbation functions as Fourier series, we 
obtain DNA rings' stability condition in a general case. By reviewing the relationship between the 
Kirchhoff model and the worm-like road chain (WLRC) model, we insert a spontaneous curvature 
term which can partly reflect the twist angle's contribution to free energy in the WLRC model 
and name this extended model the EWLRC model. By choosing suitable spontaneous curvature, 
stability analysis in this model provide us with some useful results which are consistent with the 
experimental observations. 
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I. INTRODUCTION 



An important biological function for DNA loops is that it has abundant deformations 
under different conditions. This deformations are mostly due to the interaction between 
;he DNA chain and the environment. For instance, protein operation ij, ion concentration 
2, 3 1 and temperature change all can induce deformations. Since it has been found that 
DNA shapes and deformations can be investigated by the 
work have been done based on this theory 

Baa 



experimental techniques, such as optical tweezers 
other techniques 19l . 
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Kirchhoff elastic theory, much 
, Q, 3 • However, modern 
15|, micromanipulation [l6. [itI. [isl 



3, 



11 



and 



1 

indicate that DNA has abundant mechanic characteristics 
which cannot be explained by classical Kirchhoff theory. For example, under different pulling 
force, DNA has different stretching ability 15j. Thus, several different models have been 



suggested to describe DNA chain in different conditions. For instance, the wormlike chain 



(WLC) model 



22 



model 
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was used to describe DNA under a small external force. And the WLRC 
25l | is appropriate to describe DNA with its double-helix structure under 
a moderate force. Further, following the increase of external force, new energy terms are 
introduced 26||, and the results are highly consistent with the experimental observations. 



Based on the Kirchhoff theory, the free energy density of duplex DNA chain can be 
written as 



(1) 



where Ki = — X2 ■ {d^K^/ds), K2 = ^i- (d^K^/ds), K- 
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^2 



((ixi/(is) and {xi,X2,X3} denotes 
27|, and s is the arc length of DNA 



the basis of the local coordinates on DNA chain 
centreline, oi, 02 and 03 are constants. Ki and K2 are spontaneous curvature and is 
spontaneous tension. Let = (j){s) be the angle between xi and the main normal n, we have 



Ki = K cos 0, K2 = K sin 0, = r + 



(2) 



where K = K{s) and r = r(s) are the curvature and torsion of the DNA centreline. An 
overdot denotes a differential with respect to s. Then the energy density in ([T]) can be 



2 



changed into 



J-' = J-'{K,T,(f),^) = -fC^(ai cos^ + a2 sin^ 0) 

—K{aiKi cos (j) + a2-S'2 sin > 



and the shape equations are shown in [27|]. Considering the energy density in 



+-a3(r + — Ks)^ + constant, (3) 



contains 



thus, any T{K, r) model, such as the Helfrich model [28|, the WLRC model [23|, l2J, [25 1 



and the WLC model 22|, cannot be equal to the classical Kirchhoff model, because they 
cannot give any information about the twist angle 0. But specially choose a\ = a2 = A and 
as = C, and let = c = constant, the energy density in is reduced to 

A C - 

T = -K^ - toAK + — (r + c - K^f + constant, (4) 

where to = Ki cos0 + K2 sin0. If to = 0, we get the WLRC model. But in this paper we 
choose to = constant and = c = constant, and name this model the extended worm-like 
road chain (EWLRC) model. We think the EWLRC model with to = constant is better than 
the WLRC model to a certain extent because the arbitrary constant to can partly reflect 
that the twist angle and spontaneous curvatures Ki and K2 influence DNA shapes. Then 
we choose energy density as the form 

T=^iK-tor + ^ir-coor, (5) 

where uq = K3 — c and we also name it the spontaneous tension. Clearly, the EWLRC 
model is different from the model discussed in [29|] by choosing ai = 02, K2 = and K3 = 
in©. 

By studying the second variation of DNA's free energy, Zhao et al. [28] found that the 
chiral term in energy function could induce abundant shape transformations of DNA rings 



in the Hel: 
In 
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"rich model, which is in good agreement with the experimental observations in 



291], Zhou et al. showed that the spontaneous curvature plays a significant role 
in the stability of kinked DNA. By choosing Euler angles as variable for the free energy 
function, Panyukov and Rabin |30] generally investigated the effects of thermal fluctuations 
on elastic rings and pointed out that spontaneous curvature is important in affecting the 
spatial configurations of the ring. In this paper, we discuss the second variation of the 
J-'{K, T, K, f) model and give stability condition of ring solution in Sec. [TTl In Sec. IIIH 



we take the EWLRC model as an example and obtain some useful results which are very 
consistent with experimental observations. Finally, a short discussion is given in Sec. IIVI 



II. SECOND VARIATIONS OF FREE ENERGY AND STABILITY CONDITION 
OF DNA RINGS 



We choose the free energy for the closed duplex DNA with the form 3l|, |32| 



F = y J^{K, T, k, f)ds + \jds, (6) 

where JF is an arbitrary function and ds is the element of the arclength s, and A is the 
Lagrange multiplier. The central axis R(s) of duplex DNA chain under small perturbations 
can be written as 

R'(s) =R(s)+^(s)n + y;(s)/3, (7) 

where il'is) and ^{s) are two small smooth functions, and /3 and n are the binormal vector 
and the main normal vector, respectively. Under the perturbations, the general shape equa- 

32| by studying 5^^'> F = 0. For the shape equations, a ring 



tions of DNA are attained in 



solution with radii R = 1/K and r = induces 

A = - ^0. (8) 

Here we define J-'i = ^2 = §75-^3 = *^'' ~ ^ -'^ ^^^^^ ^^^^ 

K = 1/R and r = 0, such as J-'i = J-'i\{k=i/r,t=o)- 

For different DNA models, they have different energy density functions : J^{K,t,K ^f). 
For an arbitrary function JF, we have 

^(1)^ = + 5('Vj-2 + 5^^^kT^ + 5(1)7-^4, (9) 

+ ^(5«K)^^n + \{6'^'^rfT,, + \{6^'^kfT,, + \{6'^'^ffT,,. (10) 
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The first and second variations of total energy are 



g 2 ds. 



(11) 
(12) 



To obtain S^^^F and 5(2)^, tliose variations: d'^^'^K, S'^^^K, S'^^^K, 5^^h, S'^^h, S'^^h 

and are indispensable. We give a method in Appendix A to obtain those useful terms. 

For a planar ring with K = 1/R and r = 0, we expand the two arbitrary functions ip and 
(p as the following Fourier series 

oo oo 

= Co + Cn cos{nKs) + dn sm(nKs), 

n=l n=l 

oo oo 

Co + e„ cos(?2-ft's) + /i„ sin(nKs), 



n=l 



n=l 



and after a lengthy calculation, we get the second variation, it reads 

oo 



71=1 



(13) 



with 



= 




-3 2 ttO 

115 






Pn = 




'3(^2 - 1 


)^(n^i?-^3°3 + -^i)5 




Qn = 




-V(n2 - 


- 1) [/?-i(n2 - l)(i?-V^°4 + J-°2) 


+ -^?]5 


rrin = 


2 


-2/2 

n[n — 


l)[i?-i(n2-l)(i?-V^3°4 + J-°2) 


-■^2°] 5 




2 


-4 2/2 

n [n - 


-l)'(-^2°3--^?4)- 





(14) 



It is easy to find that n = 1 is corresponding to a trivial translation. In the Eq. (fT3l) . if we 
don't know whether p„ and g„ are positive or not, it seems difficult to discuss > 0. 

Specially, if p„ and are nonnegative, such as the Helfrich model 
> yields 



stability condition 



(15) 



PnQn - ml - sl> 0. 



For the Helfrich model and any F{K, r) model, we have jFgg = JF^^ = 0, then condition (fT5|) 
is reduced to Pn9n — > 0- However, if p„ and g„ are nonnegative, only 2mn{cnhn — dnCn) 



and 2sn{cnen + dnhn) in (fT3|) have the possibihty to be negative. So any stable deformations 
need 

< or c„e„ + dnhn < 0. (16) 

It indicates that any deformations in this case are nonplanar, because we cannot choose 
= hn = 0. Then, tangential and normal deformations are coupled each other and they 
must occur simultaneously. Only when p„ < in (fT3|l can induce planar deformations 

(c„ 7^ 0, (i„ 7^ and e„ = /i^ = 0). 



III. THE EWLRC MODEL 



For the EWLRC model in ([5]) we assume A > and C > 0, then the stability condition 
( IT5|) is reduced to 

A = A{n^ -1)[A + C{n^ -1)] 

-A^{n^ - l)toR - C^uj^R^ > 0. (17) 

By solving the above inequality, one gets the stable range of R. Here, we only show the 
upper boundary as an example 



R< ' 



- iy{tl + 4W^)+4W\n^ 



+ (18) 

where W = C/A and we let c^o > 0,to > 0. For DNA rings with fixed R, such as 168- 
bp DNA with i? ^ 9.1 nm, inequality f|T8l) implies that the lower states will be instable 
and rings will change into kinked shapes if to decreases. Let to be negatively related to the 
concentration of ion, this conclusion is consistent with the experimental results that complex 
deformations will emerge following the increase of ion concentration However, we have 
reverse conclusion for ujq that it needs to be positively related to the concentration of ion. 
Specially, if c^o = 0, we need that to is positively related to the concentration of ion (we will 
see it in the late text). 
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32| 



Considering the helix solution K = ^:j^^2 ; t = j.'^^h'^ ; the shape equations in 
are reduced to (A = 0) 

A[{rl + h^)to - ro][rl{l + roto) + h\roto - 2)] 

+Cro[{rl + h^)iUo - h][{rl + h'')ojo + 3/i] = 0. (19) 



TABLE I: The parameters for helix of A-, B-, and Z-DNA. Here tq, ro//i and come from 33|, 
to is obtained by solving Eq. ([19]). 



DNA tq (nm) r^/h ujq (nm ^) tg [p-'ca. ^) 



A-DNA 1.3 3.32 2.55 0.69 or -0.56 

B-DNA 1.0 1.89 1.89 0.77 or -0.33 

Z-DNA 0.9 1.24 1.38 0.66 or 0.22 



Choosing A = 50 nm, C = 0.15 nm and the values of rQ,h,uJo in [33[, we obtain the 
corresponding shown in Table [H Specifically, letting h = 0,rQ = R, we get the equation 
for a ring 

A{1 - Rhl) - CR^ul = 0. (20) 

Solving the above equation, we get in Table HI] for 168-bp DNA and 126-bp DNA {R ^ 6.8 
nm). However, these constants in Tables [I] and [IT] only give us some coarse values, because 
the experimental observations require that uq and to should relate to ion concentration. 
Actually, using these constants is difficult to explain the phenomena that 126-bp circles in 1 
mM Zn^^ mostly are not kinked but 168-bp circles mostly are kinked in the same condition 

a. 



Specially choosing W = 0.154 and Uq = suggested by Thamwattana et al |33|, condition 
( [T8]) is reduced to 

R<^- ^ —. (21) 

^0 

Clearly, we need that to is positively related to the concentration of ion to explain that hight 
state deformation occurs following the increase of ion concentration jsl. Moreover, with the 
same W and to, Eq. ( [21]) indicates that a big DNA ring will have higher state deformation 
than a small one has. Supposing two rings with radius R\ and R2 {R\ < R2) are in the same 
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TABLE II: The spontaneous curvature to (n™ ^) obtained by solving Eq. (j20|) of A-, B-, and 
Z-DNA. 



DNA A-DNA B-DNA Z-DNA 



168-bp - 0.037 0.080 

126-bp 0.046 0.104 0.126 




FIG. 1: Stable ranges of DNA rings with W = 0.154, cjq = and different to (nm~^). From a to 
b, elliptic shapes will emerge for 126-bp DNA; from c to d, there are trigonal shapes for 168-bp 
DNA, which is consistent with the value that the number of kinks per circle is 3.83 Q]. 

condition (with the same W and to), letting T{n) = {n^ — 1)W + 1, we have 

^ < <o < (22) 

where, ni and n2 {ni < ^2) lie on the deformation level of the two rings with radius Ri 
and R2, respectively. Han et al |3l] find that 126-bp circles (-Ri = 6.8nm) in 1 mM Zn2+ 
mostly are not kinked and have somewhat elliptical shapes [rti < 3) and 168-bp circles 
{R2 = d.lnm) in 1 mM Zn^+ mostly are kinked, then we have 

Tirii) T(3) 

<to< (23) 
9.1 ° 6.8 ^ ^ 
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We find that only n2 = 3 can ensure to is existent. So trigonal deformation will emerge for 
168-bp DNA, which is consistent with the experimental result that the number of kinks is 
3.83 per circler for 168-bp DNA in [3']. Then, we get 0.245 nm~^ < < 0.328 nm~^. Fig. 
[T] shows the stable range of two kinds of DNA rings with W = 0.154 and cjq = in fl?Il) . 
from which we can see that trigonal deformation will arise for 168-bp DNA rings and elliptic 
deformation will emerge for 126-bp DNA rings. 
In f l22l) the existence of to induces 

W[inl - l)/?i - K - 1)^2] <R2- Ri- (24) 

Choosing Ri = 6.8 nm, R2 = 9.1 nm and ni = 3, we find that the above condition can 
always be satisfied with any W (note we assume W > 0) when n2 < 3, because the left 
side of fl2^ is negative. This implies that 168-bp DNA rings with any W will change into 
trigonal shapes when 126-bp DNA rings change into elliptic shapes in the same condition. 



IV. CONCLUSIONS 



We have generally investigated the second variation of DNA rings and shown the stability 
condition. However, only and if only > and gn > condition fll5l) is valid. The 
EWLRC model and the Helfrich model are satisfied with this condition, and deformations 
in those models are must nonplanar. If c^q = 0, the stability condition fl2Tl) needs that to 
is positively related to the ion concentration. Considering tg = -^1 cos0 -|- A'2sin0 and we 



choose = constant, which is consistent the conclusion in [29|] that the intercalation of 
Zn^+ takes place mainly at the exposed side of the DNA ring hence causing an increase in 
its spontaneous curvature. In Han et al. find that there are somewhat elliptic shapes 
for 126-bp DNA and that the number of kinks is 3.83 per circler for 168-bp DNA in 1 mM 
Zn2+, by choosing W = 0.154, c^o = and 0.245 nm'^ < to < 0.328 nm'^ in the EWLRC 
model, our results are very consistent with this phenomenon. 
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Appendix A 

In this appendix, we present a method to obtain S'^^^K, S^'^^K, S'^^^K, S^'^^K, ^^^V, 5^^V, 
6^^^f and For simplicity, we only give the method and key steps. The variation of 

R(s) is 

5R = R' - R = + (p(3. (25) 

Let g — R,j; -R^a; (R^^, — x is the variable in orthogonal coordinates), we have ds = s/gdx 
and 

5g'^ = -g^KiP + ^g^ [^^ ^ ^^2 ^ ^2y ^ ^2 
+2(V'0 - Mr] + 0(3), 
Sg^ = g^K^ - ^g^ [ij^ + iij^ + + 

-2V^2^2 ^ 2{iI:lp - V^(^)r] + 0(3). (26) 

Where 0(3) means the third and higher orders of ip and Lp. For an arbitrary function 
V — V(s), there is 

5Y={1 + ghg^){6y),, + g-^VSg^. (27) 
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Using above equation, we have 



5R = il + g-^5g'^){5R),s + g-^il5g^, 
SK = {1 + ghg^){6K),s + g^'RSg^. 



Then, we get 



Consequently, we attain 



6K 



6K^ = 2R 6K + 6K 6R. 



[K"^ - T^)i) -Tip - 2T(p + -11) + K 



2„/.2 



+ Klp{')p — Tip) 



Now, we have 



+ -K-\p + Ti) + 2t^ - rV)' + 0(3). 



6K-^ = -K-^6K + {5Kf + ■ ■ 



The variation of n is 



The variation of ri can be obtained by 



(5ri = {1 + g'^Sg ^ ){Sn)^s + g^iiSg 2 



The first variation of tension is 34 1 



5(i)r = K{2tiIj + p) + K-^k{T^p - fip - 2t^ - (p) 



Considering n ■ n = K + r , we have 



5(2)r = It-' 
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Inserting (129|) . (l32l) and (l33i) into (IMI) . after a lengthy calculation, we get 



+T''^/j(p) +K ^ ~ K^i){ip + tiIj) 



+2r^('?/'^ - (f^ + -ipip - (p^p - 2fipLp) 

+f(3v3v9 - 3?/;^ + v^^^V - V'^-^V) 
+2r((^2 _ ^2 ^ <^(3)0 _ ^(3)^) 

+r^f(?/'^ - ip^) + 8T^T{ijtp - ipip) 
+2ffip(p + AfT{ipLp + t/;!/? + 3'?/'0) 
+f(</3y9 — ipip) + 2f T{ip(f + ■i/'f/)) 
-■?/;^^)(^ - t/^v?*^^^ - r'^('?/;v9 + 
+2K'^k{T'^ilj + Tip + 2Tip - -ii)) 

X (r^(y9 - fV^ - 2r^ - v3). (35) 

Considering (p7|) . we have 

<5ii: = (i + (7^5^7^)(5i^),. + ^7^ii:5^?^, (36) 

5f = (l + ^55^^)(5r),, + ^5f5^^. (37) 

Here, for simplicity, we only give 6K and 6f for a planar ring with K = 1/R and r = 0, 
they are 

= + (38) 

^(^^ir = 3R~^iljip + R-\3iIj-^^^^ + 3^-^ - 2^(p) + R(p^^'^\ (39) 

5«f = R-^^ + R^^^\ (40) 
5(2)f = i?-2(3^y3 + 2^^^) - i?'(V^<^^^^ + (^V-^'^ + 2V'(^V^^^) 

+4?/^y3 + 4V'<^^^^ + 3^/'(^(^) + 2^(^V- (41) 
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